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ABSTRACT 

I review recent progress in understanding inclusive quarkonium pro- 
duction in hadron collisions. The first part focuses on non-relativistic 
QCD as an effective theory. I discuss its differences from and similar- 
ities with effective theories describing bound states of a single heavy 
quark, as far as matching calculations beyond tree-level and power 
counting are concerned. The second part summarizes predictions for 
charmonium and bottomonium production at collider and fixed tar- 
get experiments and their comparison with data. The emphasis here 
is on novel signatures due to color octet production, polarization of 
quarkonia and the Xi/x% ratio in fixed target collisions. 
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1 Introduction 



The discovery of charmonium bound statesd opened the world of heavy flavours, 
its wonderful variety and complexity. Since then, the focus has shifted to decays 
of bound states of a single heavy quark, their implications for CP violation and, 
perhaps, 'New Physics'. Quarkonia, due to their leptonic decay signature, on the 
other hand have become important tagging modes in hadron collisions and may 
become so for the quark-gluon plasma phase. For a theorist, therefore, the interest 
in quarkonia is mostly of intrinsic nature. The experimental data are there, and 
we are challenged to explain them. 

Quarkonia are the 'atoms' of the strong force. If their Rydberg energy were 
larger than A, the dynamical low-energy scale of the strong interaction, quarkonia 
would be weakly coupled bound states of a heavy quark and antiquark, tractable 
in perturbation theory just as positronium is in electrodynamics. Neither charm 
nor bottom quarks are heavy enough to satisfy this requirement. And the top 
quark decays so rapidly that toponium has barely time to form. The binding 
of charmonia and bottomonia must be described non-perturbatively. Once the 
simplicity of a Coulombic bound state is foregone, having two heavy quarks to bind 
rather than a heavy and a light quark adds complications. Apart from the mass 
scale rriQ, set by the mass of the heavy quark, a quarkonium bound state involves 
three essential 'small' (compared to tuq) scales: itiqv, the typical three-momentum 
of the constituents in the quarkonium rest frame or the inverse quarkonium size; 
m,QV 2 , the scale of binding energies and A. In a heavy-light meson, there is no 
other dimensionless ratio of scales besides A/tuq. Although none of the non- 
perturbative properties of quarkonia are calculated in the approach described in 
the sequel, the multitude of low-energy scales entails more complicated power 
counting rules than those applicable to heavy-light mesons. One may envisage 
different power counting schemes, depending on the relation of the low energy 
scales. 

The success of non-relativistic potential models in describing static properties 
of the charmonium and bottomonium family suggests that these states are in- 
deed non-relativistic and that v 2 could be used as a parameter for systematically 
expanding about the non-relativistic limit. With v 2 being small, quarkonium pro- 
duction involves two different time scales: the scale I/tuq on which a QQ pair 
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is produced^ and the scale l/(rriQV 2 ) on which the heavy quark pair binds into 
a quarkonium. Provided that the two stages of the production process can be 
separated and assuming that perturbation theory is valid at the scale tuq, the 
heavy quark production part could be computed perturbatively; anything related 
to quarkonium formation could then be factorized into quarkonium-specific, but 
production process- independent, non-perturbative parameters. 

A systematic realization of these ideas has been developed by Bodwin, Braaten 
and Lepage.i It is based on an effective field theory, called non-relativistic QCD 
(NRQCD), combined with the methods of perturbative factorization and provides 
us with a tool to calculate inclusive quarkonium production cross sections as a 
double expansion in a s and v 2 , and to leading order in A/itiq in production 
processes with (light) hadrons in the initial state. This development is summarized 
in Sect. 2. 

For phenomenology the most important insightll following from the NRQCD 
description of quarkonium production is that relativistic effects are very large in 
the production of 3 S\ states. Because v 2 is not very small, the radiation of gluons 
at late times in the production process, when the QQ pair has already expanded 
to the quarkonium size, turns out to be favored as compared to the radiation from 
an almost point-like QQ pair. As a consequence, the QQ pair can remain in a 
color octet state at distances of order 1/rriQ, a possibility that is ignored in the 
earlier color singlet model. The importance of these color octet contributions is 
supported by the large direct J / if) and if)' cross sections observed in pp collisions 
at the Tevatron.! 

Subsequent to this initial success, almost all quarkonium production processes 

have been reconsidered in the light of NRQCD. The coverage of all production 

processes is beyond the scope of this presentation and I am restricting myself to 

an overview of the current status of inclusive quarkonium production at colliders 

and fixed target in Sects. 3 and 4, respectively. A discussion of other interesting 

production processes, such as photoproduction, e + e~ annihilation, Z° or B decay, 

together with a (by now incomplete) list of references can be found in Refs.iHi. 

There is also an increasing interest in polarization phenomena, as quarkonium 

polarization provides a calculable in NRQCD, and sometimes striking signature 

*The scale l/m<g need not appear if the quarkonium is produced through a weak interaction. 
Compare B — > J/ ipX mediated by a b — > ccs transition with B c — > J/ipX mediated by a b — ► cud 
transition. In the latter case, the cc pair is produced at distances of order of the B c radius. 
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of color octet production. I believe that although the above choice of production 
processes is selective, it covers, on the one hand, the most dramatic new pro- 
duction mechanisms, and on the other hand illustrates the difficulties connected 
with a quantitative confirmation of the universality of long-distance parameters, 
as assumed in NRQCD. Especially for charmonium production, quarkonium bind- 
ing effects that are at the center of tests of NRQCD factorization can rarely be 
exhibited in isolation from other QCD effects that reside in the short- distance 
part or are neglected, such as small-rc and soft gluon effects (colliders, photopro- 
duction), BFKL-type situations (photoproduction, z — > z max ~ 1), higher- twist 
effects (fixed target, photoproduction), and parton-hadron duality (B decay). Ev- 
erything taken together makes for an intricate combination of QCD phenomena. 
Beyond any doubt, NRQCD is the correct theory for quarkonium systems in the 
heavy quark limit. Whether the charm quark mass is large enough to justify an 
expansion around this limit, will be decided by confronting predictions with exper- 
iments. As of now, we are only beginning to assess theoretical uncertainties and 
to sort out those observables that eventually will stand as solid tests of NRQCD. 

2 Effective theory and factorization in quarko- 
nium production 

2.1 Non-Relativistic QCD 

Let me assume that indeed mq ^> mQV,mQV 2 , A. In this situation of well- 
separated scales, it is intuitive that a quarkonium production cross section should 
factorize into the production cross section of a QQ pair times the probability that 
this QQ pair forms a specific quarkonium state. 

To make the scale separation explicit, one may think of integrating out all 
high momentum modes in the path integral down to a factorization scale /i in 
between mq and mqv. The result would be an effective action, which is highly 
non-local on distances 1/mq. But since the low momentum modes that dominate a 
quarkonium bound state can not resolve such small distances, this effective action 
could be expanded in an infinite series of local interactions. This expansion would 
realize an expansion in v 2 ; to all orders, it would be exactly equivalent to QCD. 
In practice, this procedure can be carried out only in very simplified situations. 
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However, as long as the strong coupling a s is small enough at the scale fi, the 
effective Lagrangian can be constructed perturbatively to a specified accuracy. 

First, one identifies the low-energy degrees of freedom. In the non-relativistic 
limit, intermediate states (in the sense of time-ordered perturbation theory) con- 
taining heavy quark pairs are suppressed and integrated out. The heavy quark 
field Q splits into a two-spinor quark field ip an d a two-spinor antiquark field x- 
The effective Lagrangian a priori consists of the most general Lagrangian, includ- 
ing non-renormalizable operators, consistent with the symmetries of QCD [J]. The 
coefficients of the operators are then tuned to reproduce QCD by comparing on- 
shell Green functions computed in QCD with those computed with the NRQCD 
Lagrangian. The desired accuracy and the values of a s (n) and v 2 determine which 
operators have to be kept in NRQCD and to what loop-order the comparison has 
to be done. Note that since the 'matching' is carried out at a scale much above 
the bound state scales, only scattering diagrams have to be computed. By con- 
struction, the NRQCD Lagrangian already has the same infrared behaviour as 
QCD. This implies in particular that the result of matching is independent of how 
the small scales ttiqV, rriQV 2 and A are related. However, their relation does have 
some consequences for which operators have to kept in NRQCD to achieve the 
desired accuracy. 

Non-relativistic effective theory has first been introduced by Caswell and Lep- 
age! as a tool to manage bound state calculations in QED. This application 
is particularly transparent, as both the short- distance matching and the long- 
distance contributions can be calculated perturbatively. The main advantage of 
non-relativistic QED compared to the Bethe-Salpeter approach is that the physics 
above the scale m e is encoded once and forever in the effective Lagrangian. Be- 
cause NRQED still contains two scales, m e a and m e a 2 , the scale separation is still 
not complete, but separating m e already entails great simplifications. See Ref.i 
for a state-of-the-art calculation in NRQED. 

The NRQCD Lagrangian takes the form 

AmQCD = £2 + £4 + £glue + (1) 

tThe NRQCD Lagrangian is usually written in the heavy quark rest frame and is therefore 
constrained only by rotational symmetry. Of course, the result of any calculation is Lorentz 
invariant. The 'hidden' full Lorentz symmetry constrains some of the coefficient functions in the 
Lagrangian. 
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The light quark and gluon part of the QCD Lagrangian remains unaltered and 
is not indicated. The contribution to the effective Lagrangian that involves two 
heavy quark fields can be obtained at tree-level from a Foldy-Wouthuysen-Tani 
transformation, generalized to the non-abelian case (see e.g. Ref.0): 



£2 = ^ 



D 2 

iDn 



2m Q 



8m% 2m 



+ -\^{D -gE-gE- D)ip + -\^{iD x gE - gE x iD)i/> (2) 
+ . • • + charge-conjugated for the antiquark. 

The leading term in square brackets describes a non-relativistic Schrodinger field 
theory. To reproduce the on-shell Green functions in QCD at order v 2 , the four 
subsequent terms have to be included. Radiative corrections due to hard gluons 
shift their coefficients away from their tree level values q = 1. The ellipses stand 
for higher order terms in v 2 . To reproduce Green functions with 2n external heavy 
quark fields, NRQCD must contain local operators with 2n quark fields. In the 
following, only operators with four quark fields are of interest. The generic form 
of £4 is 

A = E^WXx^) (3) 

i m Q 

+ four quark scattering operators, 

where Ki is a matrix in spin and colour indices and may also contain factors of 
spatial derivatives D/tuq in case of higher-dimension operators. The coefficients 
di of the annihilation operators are complex. Their imaginary parts describe the 
annihilation decay of quarkonium states.! Finally, integrating out heavy quark 
loops leads to higher- dimension operators in gluon fields, C g i ue , the non-abelian 
analogue of the Euler-Heisenberg effective Lagrangian. 

The NRQCD Lagrangian C 2 + £ g iue coincides with the Lagrangian of heavy 
quark effective theory (HQET).0 Nevertheless, the two effective theories are dif- 
ferent, because their power counting schemes are different. Because A is the 
only low-energy scale in heavy-light mesons, the importance of operators in the 
HQET Lagrangian is ordered strictly by dimension. Consequently, the kinetic 
energy operator ip^ D 2 / (2mQ)ip is suppressed by A/rnq and the leading effective 
Lagrangian ijriDQijj describes a static quark. In a non-relativistic bound state we 
expect D ~ itlqv, but Dq ~ rriQV 2 and the kinetic term can not be neglected. A 
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Figure 1: Coulomb correction to point-like ip'x creation. 



more compelling argument arises, if one begins to compute Green functions in the 
quark-antiquark sector of NRQCD. In the following I will implicitly assume the 
Coulomb gauge, which makes the physics of non-relativistic bound states most 
transparent. The amplitude for a QQ pair created in a point and then interact- 
ing through exchange of a Coulomb gluon (see Fig. [I]), the '00 '-component of the 
gluon propagator in this gauge, is given by 



%C F g\ 



d A k 1 1 1 

(4) 



(2vr) 4 p 



Po + k ~^ + ie 



On-shell po = p 2 / (2mg). The integral is done most easily by closing the contour in 
the complex fc -plane and picking up the residues of the enclosed poles. The poles 
are located at ko = po — (p + k) 2 / (2tjiq) + ie and k Q = —p + (p + k) 2 / (2itiq) — ie, 
one on each side of the real axis. Then 

T n 2 f d3k m Q C F 7ca s m Q . 

1 = C F g„ / . .„ 7=r- ^ = r^r + imaginary. (5) 

^ S J (2tt) 3 k 2 (k 2 + 2p- k-ie) 4 \p\ 

The imaginary part is divergent and related to the Coulomb phase; the real part 
exhibits the well-known Coulomb divergence close to threshold. 

In the static approximation, the 'kinetic term' in the propagators in @ would 
be dropped and the on-shell condition reads po = 0. In this limit the integration 
contour becomes pinched between the two poles and the integral becomes ill- 
defined. Thus, the kinetic term must be kept in the propagator to regulate the 
pinch-singularity. If one of the quarks were light, as in physical processes to which 
HQET applies, the static limit can be taken. Although the static propagator pole 
lies on the real axis, the contour can be deformed away from the pole. 
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From (H) it follows that the effective coupling for the exchange of Coulomb 
gluons is a s rriQ/\p\. At small momenta \p\ ~ mQa s , Coulomb exchange can not 
be treated as a perturbation, even at weak coupling a s . The result of resumming 
Coulomb gluons is of course well-known and leads to Coulomb bound state poles 
in the quark Green function. On the other hand, HQET contains only strongly- 
coupled bound states. Note that, for the purpose of matching, it is sufficient to use 
free NRQCD propagators. Being of infrared origin, the Coulomb-enhanced terms 
cancel in the matching coefficients at every order. The integral in @ contains 
only one scale \p\ = vtiqv. Thus, its finite real part is dominated by gluons with 
momenta vtlqv 

Because the leading order Lagrangian contains rriQ in NRQCD, NRQCD does 
not lead to flavour symmetry. The non-perturbative properties of charmonia and 
bottomonia remain unrelated even in the non-relativistic limit. On the other 
hand, the NRQCD Lagrangian exhibits heavy-quark spin symmetry at leading 
order in v 2 , which turns out to be useful to reduce the number of non-perturbative 
parameters that describe quarkonium binding. 

I started out with quarkonium production, but the effective Lagrangian above 
does not yet allow quarkonia to be produced. One could not have expected that, 
because the short- distance part of a production process depends on the initial 
particles (hadrons, photons, Z bosons ...) that initiate it. But for every production 
process, the separation of short and long distances can be performed in the same 
way as for the effective Lagrangian above. I will elaborate on this in the following 
subsection. 

At this point, let me pause for a discussion of the scales rriQ, rriqv and rrigv 2 in 
charmonium and bottomonium systems. They are shown in Tab. [I], together with 
those for positronium for comparison. The values of v 2 are based on potential 
models, which describe the spectrum of quarkonia reasonably accurately. For 
positronium v ~ a. In a quantum field theoretic context, v 2 could be defined 
as the expectation value of a derivative operator that scales like v 2 according 
to the scaling rules to be described later. In general, v 2 will be a complicated 
function of itlq and A. In the limit vtlq — > oo, the binding becomes Coulombic 
and v ~ 1/ ln(mQ/A). Neither charmonia nor bottomonia are Coulombic, because 
the energy scale vtiqv 2 is of order A for both quarkonium families. The same 

•f However, they can not be thought of as on-shcll intermediate states in time-ordered perturba- 
tion theory and in this sense they are not 'dynamical'. 
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cc 


bb 


1 

e^e 




1.5 GeV 


5 GeV 


0.5 MeV 


rriQV 


i ou Me V 


1.4 Lre V 


6. 1 Ke V 


m Q v 2 


400 MeV 


400 MeV 


25 eV 


V 2 


0.25 


0.08 


5 ■ 10~ 5 


a(m Q )/7i 


0.1 


0.07 


2 ■ 10~ 3 



Table 1: Scales in onium systems and the expansion parameters of the non- 
relativistic approximation. 

conclusion can be obtained from the observation that with r ~ l/(rriQv), the 
contribution of the linear term in the Cornell potential^ 



V(r) 



(6) 



is non-negligible with respect to the Coulomb term. (The string tension is 
a w 430 MeV and C F = 4/3.) Let me stress that NRQCD does not rely in 
any way on a Coulombic system. It is enough that vhq is large compared to all 
bound state scales. The values in Tab. [1] should be understood for the ground 
state of each family. Excited states are more non- relativist ic and at the same 
time even less Coulombic. It is possible that the hierarchy of vtiqV and rriQV 2 with 
respect to A changes as one considers higher excited quarkonium states. Different 
power counting rules would then apply to different members of the onium family. 
Comparing the different onium systems, note that relativistic corrections (gov- 
erned by the parameter v 2 ) are exceedingly small for positronium, comparable 
to radiative corrections (governed by the parameter a s (rriQ)/7r) for bottomonium 
and definitively large for charmonium. 



2.2 Factorization and matching 

The factorization of a quarkonium production process, suggested in Ref.i, begins 
with the observation that the creation of the QQ pair requires an energy larger 
than 2rriQ and therefore some of the propagators of the diagram are off-shell by 
at least tuq, much larger than the typical off-shellness in a quarkonium bound 
state. These propagators can be 'contracted to a point', and the remainder of 
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Figure 2: Trivial example of factorization. 

the diagram, sensitive to distances larger than 1/mg, is then associated with an 
operator matrix element in NRQCD. This is true, provided the production process 
is inclusive, 

A + B ^ij(P,\)+X, (7) 

where X denotes light hadrons and A the quarkonium polarization state. If it 
is not inclusive, the process is sensitive to the details of the hadronic final state 
and its long-distance contributions can not be absorbed in quarkonium-specific 
objects alone. 

The short-distance coefficients can be computed by replacing the quarkonium 
in the final state by a perturbative QQ state at small relative momentum, because 
the short- distance coefficient by construction does not depend on those effects that 
bind such a quark state into the quarkonium. Hence one can (and should) use 
on-shell quark states and, instead of (^), one considers 

A + B —> Q(p)Q(p) + gluons, light quarks, (8) 

where 

P = (\/ m Q + q?, Qi) P = (V m Q + 9a> ( 9 ) 

in the frame where the QQ pair is almost at rest. (We may think of this frame 
as the quarkonium rest frame.) The amplitude squared is then expanded in the 
small quantities <fi, g*2 and the external heavy quark spinors are expressed in terms 
of two-component spinors £ and 77. For a quarkonium not at rest, one applies a 
Lorentz boost to the vectors p and p. For the case g*2 = ~ Qi explicit formulae for 
the reduction of the amplitude to rest frame two-spinors can be found in Ref.0. 
The simplest example of tree-level matching is shown in Fig. ^. Evaluating the 



9 



cut diagram, summing over all polarizations results in 
lQir 3 a 2 



5(s - 4m 2 Q ) x ^VT^eVT^ + 0{ql q 2 , ?i • ft), (10) 



27(2m Q ) : 

where s is the cms energy of the qq pair. The spinor product can be identified 
with the tree-level evaluation of the matrix element 

(Ot(%)) = £<0| X VT<V#(A) +X)^(X) + X\^a i T\\0) (11) 



where ip is again replaced by a QQ pair at small relative momentum. Eq. ((TTJ) 
displays the structure of a quarkonium production matrix element. In general, the 
QQ pair can be in various color and spin states. Evidently, in our example, the 
QQ pair must be in a color octet and spin-one state as expressed by the product 
<j l T a in flTTI). Further terms in the expansion in q\ and g*2 can be associated with 
operators similar in form, but with derivatives. It is natural to combine these 
into 'relative' and 'cms' derivatives, such that the following identifications can be 
made 



(fc + fe)*^ - ^(VV)- (12) 



Operators with one relative derivative on each fermion bilinear can be decomposed 
as 3 Pj with J = 0, 1, 2. It is important that the quantum numbers of the QQ pair 
in the operator need not coincide with the quantum numbers of the quarkonium 
ip, because they refer to the QQ state at a time r ~ 1/tuq (in the quarkonium rest 
frame), long before the quarkonium is formed. In between gluons with energies 
less than mq can be emitted (so that X becomes non-trivial) and change the color 
and spin of the QQ pair. By the nature of factorization, these low energy gluons 
have to be and are included in the definition of the non-perturbative parameters 
above. 

Operators with cms derivatives are usually neglected, because they are sup- 
pressed in v 2 according to the power counting rules discussed below(§]. These 
operators are important to resolve an ambiguity that arises at leading order in 
v 2 and which is apparent in (|10D: the phase space restrictions (s = ^itlq) are ex- 
pressed in terms of partonic variables, such as the quark mass, and do not reflect 

§For quarkonium decays such operators were introduced by Mannel and SchulerEi 
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the physical phase space for a quarkonium state with mass different from 4mg. It 
must be like this, because the phase space is part of the short-distance coefficient 
and therefore can depend only on short- distance parameters. It is also consistent, 
because — 2m,Q ~ rriQV 2 and the mass difference can be neglected at leading 
order in v 2 . This negligence is certainly unjustified, if an observable is sensitive 
to the kinematic boundaries of phase space and leads to large ambiguities even in 
fully inclusive (and pt-hitegrated) hadroproduction cross sections at high energies, 
because of the steep rise of the gluon distribution at small x. In the example of 
Fig. | qi + q 2 enters the phase-space delta-function. Expansion in q\ + q\ leads 
to a series of higher- dimension operators with increasingly singular coefficients, 
schematically written as 

Y,c n 6 {n \s - Am 2 Q ) (x^ l T a ^) (D) n (^VT a X ), (13) 

n 

where the superscript on the delta-function denotes derivatives. The resummation 
of this series leads to a (non-perturbative) distribution functionlll with support 
properties such that its convolution with the short-distance coefficient reproduces 
the physical phase space constraints^- We will meet a particular application of 
this in Sect. |3~T| . 

So far I have discussed only tree-level matching. Factorization becomes non- 
trivial, when one computes radiative corrections, for instance to the leading-order 
qq annihilation process in Fig. 

q + q^Q(p)Q(P)+g- (14) 

When the QQ pair is in a P-wave state, the cross section is infrared divergent, 
when the gluon is soft.0 (I assume that collinear singularities from emission off 
the initial quarks have already been absorbed into redefined partons.) For a long 
time, this infrared divergence has shed doubt on the perturbative calculability 
of P-wave quarkonium decay or production. Within the NRQCD approach, this 
problem finds its natural solution:0 the infrared divergence indicates that the soft 
gluon emission is sensitive to bound state scales and therefore should be factored 
in a NRQCD matrix element. Since before gluon emission, the QQ pair is in a 
color octet, spin-one state, the only candidate is (Og^Si)). Indeed, since this 

^In the context of quarkonium dec ays, t he breakdown of the NRQCD expansion near boundaries 
of phase space is discussed in Refs.ll3iIIl. 
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matrix element appeared at leading order (Fig. and fllOD), a complete match- 
ing calculation includes the a s correction to this matrix element, computed in a 
perturbative QQ state at small relative momentum. This contribution is both 
ultraviolet and infrared divergent. The ultraviolet divergence can be absorbed 
into a renormalization of (Of ( 3 <Si)), which becomes explicitly factorization scale 
dependent. The infrared divergence cancels the infrared divergence in the calcu- 
lation of the process (|l4j) . The short-distance coefficient of the P-wave matrix 
element is now infrared finite, but contains a scale-dependent ln(mg//i). The 
scale-dependence cancels with the scale-dependence of the octet matrix element 

Note that before the advent of NRQCD, the infrared logarithm lnmg/A was 
sometimes treated as an adjustable non-perturbative parameter, see for instance 
Reffi When this is done consistently in all production processes, this procedure 
is fully equivalent to taking into account the color octet contribution (Of( 3 Si)) in 
NRQCD. It is for this reason that the leading color octet contributions in NRQCD 
for P-wave production do not lead to significant differences compared to the color 
singlet mo deli at next-to-leading order. Ironically, color octet contributions turn 
out to be more important for S'-wave quarkonia, precisely because they are sup- 
pressed in v 2 and not intertwined with the leading color singlet term through 
(logarithmic) operator mixing. 

It is quite instructive to go into some details of the NRQCD part of the match- 
ing calculation. As an example, let me consider the one-loop mixing of (Of( 3 Sx)) 
into (Of ( 3 P )). To get a non-zero contribution, a transverse gluon must be cut. 
The dimensions work out correctly, because a transverse gluon couples propor- 
tional to p/irtQ. One of the diagrams is shown in Fig. |3|; it leads to an integral 

d 3 k 1 1 p ■ p' — (p ■ k)(p' ■ k)/k 2 

(15) 



(2tt) 3 2|fc| m 



I I 2mQ 



1 1 2mQ 



The integral is divergent and has to be regularized. One would like to use dimen- 
sional regularization, because it makes matching calculations particularly simple. 
Dimensional regularization is fine, as long as the effective theory contains only 
one low-energy scale, such as HQET. Then all integrals in the effective theory are 
scaleless and vanish identically. NRQCD integrals such as (p)|) are not of this 
type. 

To see how such integrals are evaluated, let me consider the following simplified 
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version of (0), both with a cut-off and dimensional regularization: 



Ac 7 1 2 

dk rrir. 



x 



k (m Q +p + k) 



2 



oo „ 2 



Jd "^ y ^(m Q+P +^- (16) 

A 

I will only consider the logarithmically divergent and finite contributions in the 
limit of small infrared cut-off A. Because NRQCD is an effective theory below 
the scale uiq, the ultraviolet cut-off must be chosen such that p < A c < itlq. 
Therefore, taking the integral J c and expanding it in A c /mq and p/rriQ, one gets 

^4-^,4-^+.... (17) 

V m Q/ X m Q 
The same result would have been obtained, if the integrand had been first ex- 
panded in p/triQ and k/mq, which are both small for k < A c . Note the loga- 
rithmic term in the cut-off A c , which would correspond to mixing of (Of ( 3 Si)) 
into (Of ( 3 P )) in the real case. The dimensionally regulated integral, in the limit 
e — > 0, evaluates to 



(m Q +p) 2 V A 

and seems to have no ultraviolet divergence. The problem here is that dimensional 
regularization does not know about the physical requirement pCA c C rriQ, that 
expresses that the factorization scale is smaller than rug, and treats the cut- 
off as if it were larger than tuq. Indeed, up to power-like cut-off dependence 
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coincides with the result for Ja c evaluated for A c 3> tuq. We can force dimensional 
regularization to treat tuq larger than all other scales by expanding the integral 
in mq before integration. The result is 

Jt = ( l -^.)(i-^) + ..., (19) 



in agreement with ( |T7D up to power-like terms in the cut-off which are always 
zero in dimensional regularization. Thus, in dimensional regularization one must 
expand the integrand before integration, while with a cut-off integrating before 
or after expansion yields the same result, provided A <C tuq. The calculation 
may be simplified even further, when the infrared divergence is also regulated 
dimensionally. After expansion of the integrand, all integrals are scaleless and 
vanish in dimensional regularization. This technique has been used in Ref.il to 
obtain the NRQCD Lagrangian up to order a s /rriQ in the single heavy-quark 
sector. 

Expanding the integrand works for matching calculations to all loops in the 
single heavy-quark sector, but fails in the quark-antiquark sector. Indeed, we know 
already from ([5]) that some matrix elements in the effective theory must be non- 
vanishing, so that the Coulomb divergence cancels in the matching. Integrands 
containing Coulomb singularities can not be expanded in mq before integration 
over k . While at one- loop it is easy to separate the Coulomb contributions 
explicitly, a general matching scheme based on dimensional regularization that 
would treat diagrams with mixed Coulomb and transverse gluon exchanges has 
not yet been devised. Since the problem does not appear in matching calculations 
with a cut-ofljjj], it is not related with the effective NRQCD Lagrangian per se. 
See, however, Ref.0 for an alternate view of the problem. 

Summarizing the discussion of this subsection, the differential quarkonium 
production cross section in a hadron-hadron collision, A + B — > ip{P, A) + X, can 
be factorized as 

da = J2 dxxdx 2 fi/AMfj/ufa) d °i+j^QQ[n]+x (Cn (A) )- (20) 

This equation is diagrammatically represented in Fig. [|, where I have replaced 
hadron B by a virtual photon for graphical simplicity. Each factor in (^) cor- 
responds to a subgraph in Fig. (| and a momentum scale that dominates this 

UnRQED calculations are naturally done in such a scheme, see e.g. Ref.i. Otherwise, the bound 
state diagrams would also have to be evaluated in dimensional regularization. 
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Figure 4: Diagrammatic representation of factorization in 7* + A —> vp + X. For 
the cross section the diagram has to be cut. 

subgraph. The amplitude squared for A — > i + remnant J is given by the parton 
distribution fi/A{%i)- The typical virtualities in this subgraph are of order A 2 . 
Parton i then participates in a hard collision H, in which a QQ pair in a certain 
state n is produced. This involves energies of order rriQ. The QQ pair (and, in 
general, additional gluons) connects to the quarkonium subgraph, represented by 
the NRQCD matrix element (O^^), which contains all soft lines that are sen- 
sitive to the bound state scale rriQV 2 . After factorizing initial state singularities 
into redefined parton distributions and performing the NRQCD matching as de- 
scribed above, the short-distance cross section da i+ j^QQ^ +x is infrared finite, but 
depends on the collinear and NRQCD factorization scale. The scale dependence 



cancels in the product (|20|) . 

Note that the factorization formula above implies complete factorization be- 
tween final and initial state; all lines connecting both run through the hard part of 
the diagram. The cancellation of soft gluons, that could connect the quarkonium 
and remnant jet part, is believed, but not yet really proven, to hold perturbatively 
for inclusive quarkonium production. This is quite intuitive, if the quarkonium is 
scattered at large transverse momentum with respect to the colliding hadrons. It 
holds also for the total quarkonium production cross section even though the bulk 
cross section is from quarkonia at small transverse momentum. This is intuitive 
only in the Coulombic limit itiqv 2 3> A, in which quarkonium formation has ter- 
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minated long before the time-scale of a second interaction with the remnant. Since 
perturbative matching is insensitive to the relation of low energy scales, factoriza- 
tion should thus hold in perturbation theory. The calculation of next-to-leading 
corrections to total P-wave production cross sections^ confirms factorization at 
the one-loop level. The mathematical statement of factorization is thus that all 
corrections to (|20D scale as some power of l/mq, or another kinematic scale larger 
than rriQ, as vtlq — ► oo. (Recall that relativistic and radiative corrections scale 
as 1/ In rriQ in this limit.) Corrections should then be suppressed as A divided 
by potentially any other scale in the process. Since rriQV 2 ~ A for charmonium 
and bottomonium, one may expect large 'higher- twist' corrections in fixed-target 
collisions, when the heavy quark-antiquark pair moves parallel with a remnant jet 
and remains in its hadronization region over a time 1/A in the quarkonium rest 
frame. Even if higher-twist effects scale only as A 2 /m 2 for charmonium, they can 
be expected to be non-negligible for total cross sections. 

2.3 Power counting (velocity scaling) 

The factorization formula (^) contains an infinite series of non-perturbative pro- 
duction matrix elements (O^^) and would be useless, if it could not be truncated 
after a finite number of terms. A first indication comes from the dimension of the 
operators, since any power of 1/tuq in the coefficient of an operator can be com- 
pensated only by one of the low-energy scales vtiqV , mqv 2 or A. However, the 
matrix elements can have additional suppressions beyond their dimension, be- 
cause of the particular structure of a non-relativistic bound state. 

The standard power counting (or 'velocity scaling') rules are due to Ref.@. For 
example, heavy quark fields scale as (mgt>) 3//2 , because / d 3 rip^ip counts heavy 
quark number and r ~ l/|p| ~ l/( m Q t 0- The Virial theorem relates the potential 
to the kinetic energy E ~ vtiqv 2 and leads to gE ~ ttiqV 3 for the electric field 
inside a quarkonium. From the equation of motion for the vector potential, gB ~ 
rnqV 4 follows for the magnetic field. As a consequence the 'dipole interaction' 
ip^(gA ■ d/rriQ)ip and 'magnetic interaction' ip^(gB ■ a/(2mg))^ in the NRQCD 
Lagrangian both scale as v 2 relative to the Coulomb interaction. 

It follows that a matrix element such as (O^ ( 3 Si)) scales as v 3 , because a QQ 
pair in a color singlet 3 <Si state overlaps with the leading Fock state wavefunction 
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of a J /if; without need of soft gluon emission^. In the following the scaling 
of all matrix elements will be considered relative to [Of ( 3 Si)}, that is, I put 
(Oi^ ( 3 Si)) ~ 1. To estimate the scaling of other matrix elements, the multipole 
suppression of gluon emission with momentum k ~ itlqv 2 <C p ~ uiqV has to be 
taken account as a consequence of the two low energy scales present in NRQCD. 
(In the context of NRQED, the separation of photons with momenta \k\ ~ m e v 
and \k\ ~ m e v 2 and the power counting for the two momentum regions has been 
analyzed in detail by Labelled) 

Consider the P-wave matrix element (0%^ ( 3 Po)) in a J /if} state. A non-zero 
overlap requires the emission of gluons into the final state. For a single gluon 
emission, one such diagram looks exactly like Fig. The QQ vertex carries a 
derivative for the P-wave operator and the transverse gluon also couples propor- 
tional to p/rriQ from the dipole interaction term above. Thus (Og ( 3 Po)) scales as 
p 4 /mQ ~ v A relative to (Of '( 3 Si)) , to which diagrams with out radiation would 
contribute. 

Consider now (O^^^Sq)), which requires a spin-flip transition. As the chro- 
momagnetic interaction is proportional to the gluon momentum, a diagram such 
as Fig. | leads to (cf. ([15])) 

Ac 



a. 



d 3 k 1 Pij (k x a)i (k x a)j 



Am 2 Q J (2tt) 3 2|fc| 



i i Pi (p+fc) 2 , ■ 



P'o+\k\-^ + ie 



A 2 

O-szpt ma g ne tic dipole contribution, (21) 



where Py = 5ij + kikj/k 2 . In the second line I have (schematically) separated 
the contribution from the region \k\ ~ vtlqv. In this region the |fc|-terms in the 
quark propagators dominate, the integrand becomes independent of the bound 
state structure and results in a pure (power-like) cut-off term. The contribution 
from | A; | ~ mqv 2 is denoted as 'magnetic dipole' and scales as 

^L\k\ 2 ~ v *\^ v \ (22) 
m Q 

where for \k\ ~ rriQV 2 the coupling a s must be counted asl. A = k/p ~ v 

corresponds to the ratio of the size of the quarkonium and the wavelength of the 

**Up to corrections in v 2 , this matrix element coincides with the wavefunction at the origin 
squared. Keeping only this contribution, NRQCD reproduces the color singlet model for S'-wave 
quarkonia. 
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(Of( 3 ^)) 


<<^( 3 Si)> 


(Ot(%)) (Ot( 3 Po)) 


NRQCD 


1 


V 4 


v 2 X 2 v A 


CEM 


1 


1 


1 v 2 



Table 2: Velocity scaling in NRQCD and the color evaporation model (CEM) 
relative to (Of (^O) for S -wave quarkonia. The 'standard' velocity counting is 
recovered for A = v. 



emitted gluon and provides the expansion parameter for the multipole expansion. 
In general, once gluons of momentum itiqv are separated, the interaction vertices 
of NRQCD can be multipole-expanded. This is true beyond perturbation theory 
and justifies the single-gluon approximation which I used for illustration. Thus 
{0^(%)) scales as v 4 , if the cut-off term in (]2l|) could be discarded. If one 
identifies A c ~ rriQV and a s at the scale rriQV with v, as suggested by a Coulombic 
limit, the cut-off term scales as0@ v 3 . Being pure cut-off, however, it should 
not be taken into account to estimate the scaling of low-energy matrix elements. 
It would be cancelled with a contribution to the coefficient function, if it were 
evaluated with the same cut-off. In dimensional regularization, the contribution 
from \k\ ~ vhqv gives a tadpole-like integral and would be set to zero naturally. 

The multipole suppression is effective as long as rriQV ^> A ~ rriQV 2 holds. If, 
on the other hand, rrigv ~ A ^> itiqv 2 , the typical momenta of (soft) gluons would 
be A and the multipole expansion would not work, as A = 1. We can keep A as 
a free parameter and conceive an intermediate case for J/ip or, in particular, if)'. 
The power counting for the most important QQ states is summarized in Tab. |2| 
and 3 for S- and P-wave quarkonia, respectively. For each QQ state, there exist v 2 
corrections due to operators with two and more derivatives on a single bilinear of 
fermion fields. Since their coefficient functions are not enhanced by fewer powers 
of a s compared to the leading operator in each channel, I will not discuss them 
further. Like operators with cms derivatives, these operators can be important in 
specific kinematic regions. 

Before closing this section, I would like to mention the color evaporation 
modelil (CEM), the only remaining potential competitor of NRQCD [tt| In the 

alluded to earlier, the color singlet model has been swallowed by NRQCD and lost its 
justification hence after. 
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(Of( s Pj)) (Of( 8 £i)> (Ot{*Pj,)) {OiQ-Px)) (Cti^i)) (Ot( 3 Dj„)) 


NRQCD 
CEM 


v 2 v 2 v 6 f 4 A 2 v 6 v 6 
v 2 1 v 2 v 2 1 f 4 



Table 3: Velocity scaling in NRQCD and the color evaporation model (CEM) 
relative to (Of ( 3 Si)) for P -wave quarkonia with total angular momentum J. The 
'standard' velocity counting is recovered for X = v. 

CEM, the NRQCD expansion on the right hand side of (|20|) becomes replaced by 
an average, 

2rriQ q 

(Ot^)^U / dMda i+J ^ mM)+x /dM, (23) 

n n J 

2rriQ 

i.e., the open heavy quark cross section is integrated over the invariant mass M 
of the QQ pair up to the open heavy flavour threshold. It is then assumed that 
the quarkonium production cross section is a universal (for each i/j) fraction of 
the sub-threshold cross section. 

Spiritually, the CEM is close to NRQCD in that it also allows a QQ pair in a 
color octet state in the hard collision to hadronize into a quarkonium. It is also 
clear that, because the NRQCD expansion arises from an expansion of the open 
QQ production amplitude at small relative QQ momentum and because rriQ q — 
rriQ <C itiq, the CEM is very similar to NRQCD as far as kinematic dependences of 
the production cross section are concerned. The difference arises in the importance 
that is assigned to the various terms that arise in the expansion close to threshold. 
In terms of NRQCD matrix elements this difference can be summarized by the 
statement that the power counting implied by the CEM would assign v 2d ~ 6 to 
any dimension d operator, independent of the color and spin state of the QQ pair 
(see Tab. [2] and |3|). The usual argument is that the emission of soft gluons in 
the hadronization of a QQ pair randomizes spin and color, so that by the time 
the quarkonium forms, any information on the state of the initial QQ pair has 
been lost. The problem with the argument is that soft gluons do not flip a heavy 
quark spin easily, a piece of information that is incorporated in NRQCD via spin 
symmetry but disregarded in the CEM. Moreover, the CEM treats a quarkonium 
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as if its size were that of a light hadron, in which case the multipole expansion is 
not valid. The CEM also does not incorporate the heavy quark limit, in which the 
binding becomes Coulombic. This does not impede the color evaporation model 
to provide a description that is phenomenologically successful in some cases, in 
particular as - in contrast to the color singlet model - the correct kinematic 
dependences are incorporated. In general, with only one parameter the CEM 
is too restrictive. It predicts are universal x/(^/V)- ra ti° fx/ fj/ipi which is not 
supported by the comparison of quarkonium production in fixed target collisions, 
photoproduction and B decay. Note that - as NRQCD - the CEM is based on a 
leading-twist approximation to the heavy quark production cross section. 



2.4 Quarkonium Polarization 

NRQCD factorization also predicts the polarization of the produced quarkonia, al- 
though, in general, at the expense of introducing further non-perturbative matrix 
elements that do not appear in the cross section summed over all if) polarization 
states ('unpolarized cross section'). Let me sketch the decomposition of matrix 
elements for the case of a S-wave, spin-triplet quarkonium.0SliJ 

After expansion of the QQ production amplitude squared in relative momen- 
tum, and after decomposing color indices into color singlet and a color octet term, 
and spin indices in a spin singlet and triplet term, the cross section can be written 
as a sum (compare with (p0|)) 

E^WRf), (24) 

n 

where the short- distance part C is still coupled to the matrix elements through 
three- vector indices that arise in the expansion in relative momentum. To decouple 
these indices, one writes down the most general decomposition of the Cartesian 
tensor (Of^) incorporating rotational invariance and the fact that the rest frame 
matrix elements can depend only on the polarization vector (in general, tensor) 
of the quarkonium. 

The constraints from spin symmetry are incorporated as follows. Because of 
spin symmetry, spin and orbital angular momentum are separately good quantum 
numbers at leading order in the velocity expansion. The angular momentum part 
of the quarkonium wave function is thus a direct product S x Li (i may be a multi- 
index for orbital angular momentum greater than 1) in spin and orbital angular 
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momentum and can be represented as 



% e a (X)a a 

l Pi ei(X) (25) 
3 Pj ^(JA|lp;la)e,(p)e Q (a)a a 

in the quarkonium rest frame, where (JA|lp; la) denotes a Clebsch-Gordon coeffi- 
cient and e(A) an angular momentum-one polarization vector. A general NRQCD 
matrix element can then be written as 

J2(0\ X ^D in . . . A n V#(A) + X) (V(A) + X\^ K 'D jn . . . D jmX \0) 
x 

= H{kS) tr(K'S) LiLj Cijh...i n n...j m , (26) 
where k is a matrix in spin and color indices and 

boSkiS mn + bi5 km 5i n + b 2 S kn 5i m (27) 

etc.. a, bi are scalar non-perturbative parameters. Some of them can be expressed 
in terms of those that appear in unpolarized production by taking contractions 
or summing over polarizations. For a >S-wave quarkonium, the orbital angular 
momentum part is trivial and its is straightforward to obtain 

£<°lxVW(A) +X)M\) + X\tfo b X \0) = \{Ot( 3 S 1 ))e a *(X)e h (X) 
x 6 

^2(0\x^ a T A 4jmX)+X)^(X)+X\^a b T a X \0) = i(Gf( 3 £i)> e a *(X)e\X) 
x 6 

£(0| X VT a Bi) ^(X)+X)(HX) + X\^a b T A {-'- B,) x|0> (28) 

= (Ot( 3 P ))5 ij e a *(X)e\X) 

^(0\x^mX)+X)(i,(X) + X\^T\\0) = \(OtCS )). 
x 6 

The matrix elements on the right hand side include implicitly a sum over po- 
larization as in the definition of Ref.i and in (|Pij). For the last line zero would 
be obtained in the spin-symmetry limit. The factor 1/3 here follows from rota- 
tional invariance. Note that the decomposition of P-wave operators in the third 
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and fourth line is not diagonal in the angular momentum basis JJ Z . As a conse- 
quence, for S-wave production, the total angular momentum of an intermediate 
quark pair in a P-wave state is not a good quantum number. In this basis, in- 
terference term between states with different total angular momentum have to 
be included to obtain a factorized form of polarized production cross sect ions .0 
These interference terms vanish when all polarizations are summed over. 

Notice that after applying spin symmetry, no new non-perturbative matrix el- 
ements had to be introduced to describe polarized ^-wave quarkonium production 
as compared to unpolarized production at this order in the velocity expansion. 
For P-wave quarkonia, this is also true at leading order in v 2 , but - in contrast 
to S'-wave states - no more for v 4 corrections. 



3 Quarkonium production at the Tevatron 

In this present section I discuss the comparison of predictions for S-wave pro- 
duction with Tevatron data,Bil which inspired to a large extent the theoretical 
development described above. In Sect. |] fixed-target data will be revisited in the 
NRQCD approach. 



3.1 Cross sections 

Let me follow the chronology of the development of theory and data. When CDF 
measured, for the first time, separately J /if) and if>' production not coming from 
B decay ('prompt'), they found much larger cross sections than expected.0 The 
expectation then was based on the color singlet process 

g + g^ccfS?]+g: (29) 

Pt 

for the direct production of .S-wave states, with additional contributions to J /if) 
from radiative decays of x& an d Xci- As indicated, this lowest order process leads 
to a very steep pt-spectrum da /dp 2 of the short- distance cross section. On the 
other hand at pt 3> IrriQ the quarkonium mass can be considered as small and the 
inclusive production cross section, like any single-particle inclusive cross section 
in QCD, should exhibit scaling: da /dp 2 ~ 1/p 4 up to logarithms. Such processes 
can be described by convoluting a (properly factorized) parton production cross 
section with a fragmentation function D^h^z), where z denotes the fraction of 
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i's momentum transferred to the hadron H. Braaten and YuanH noted that for 
quarkonia the dependence of fragmentation functions on z can be calculated in 
the NRQCD approach. Since the hadronization of the QQ pair takes place by 
emission of gluons with momenta of order rriQV 2 in the quarkonium rest frame, 
the energy fraction of the quarkonium relative to the fragmenting parton differs 
from that of the QQ pair only by an amount 5z ~ v 2 <C 1. They argued that the 
higher-order color singlet process 



g + g^ cc^j+gg + g : a]- A (30) 

Pt 

exceeds (|29| ) for p t > 7 GeV as relevant to most of the Tevatron data. This 
expectation was born out by detailed analyses.0 At this point, taking also into 
account fragmentation production of P-wave states and their subsequent decay, 
theory and data seemed to agree for J/ip production, but showed a factor of 
30 discrepancy in overall normalization for if)', a discrepancy then referred to 
as ^'-anomaly. The discrepancy in normalization indicated that an additional 
fragmentation contribution had been missed. Braaten and Fleming! suggested 
that the color octet fragmentation process 

g + g^cc[ 3 s{ 8) ]+g: a 3 s v*\, (31) 

Pt 

where a gluon fragments as g — > cc[ 3 S[ 8 ^}, could make up for the missing piece. 
The price to pay is a new parameter (Of ( 3 Si)) that had to be fitted to the data. 
Because (|3l|) has two powers of a s /ir less than fl30D, the magnitude of (Of ( 3 5i)) 
turned out to be consistent with its v 4 suppression according to the scaling given in 
Tab. 0. Subsequently, CDF was able to remove x feed-down from J/ip production. 
The same discrepancy as for ip' appeared in the J/tp cross section and could be 
explained as for ip' by color octet fragment at ion. @ 

With this explanation of the ^'-anomaly at hand, further studies focused on 
finding additional consistency checks of the color octet hypothesis beyond the con- 
sistent size of the octet matrix element. One such check derives from polarization 
and will be described in the next subsection. Others come from different charmo- 
nium production processes. Apart from fixed-target production discussed in the 
following section, I do not undertake such a comparison in this article, a compar- 
ison that would in any case be still preliminary. Let me instead continue with the 
Tevatron analysis and discuss the uncertainties in the theoretical prediction. 
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At moderate pt ~ 2m,Q two further octet production channels, which do not 
have a fragmentation interpretation at order a 3 , 

g + g^ce£si?\ 3 PP)+g: (32) 

need to be considered, too, even though they seem to be suppressed by v 4 with 
respect to the color singlet process (p9|) . (I'll come back to this point later.) These 
contributions have been calculated in RefsJllIzI and turn out to be significant for 
p t < 10 GeV. In this region (and taking m c = 1.5 GeV) the transverse momentum 
distribution is sensitive only to the combination Mf^S^^P^) defined as 



M^(^ 8) , 3 P (8) ) ee (OfCSo)) + — 2 (Ot( 3 P )). (33) 

mi 



This matrix element appears as a second fit parameter. The fits to the CDF 
dataH are shown in Figs. |5] and || These fits are based on a combined fit of J/ ip 
and ijj' data, keeping the ratios of octet matrix elements for J / ip and tp' fixed. See 
Ref.0 for further details. The values of the matrix elements with parameters as 
specified for Fig. |5] are found to be 

(Og^f^)) = 1.06- 1(T 2 GeV 3 

(O^fSo)) + ^4 (Oi^CPo)) = 4.38 ■ 10- 2 GeV 3 

(Of ( 3 Si)) = 0.44 • 10" 2 GeV 3 (34) 

(Of? So)) + ^1 (°t'( 3 Po)) = 1-80 • 10- 2 GeV 3 . 

Note that Figs. [| and || include only a 3 processes. The leading color singlet process 
at large pt, Eq. (^), is not included but would fall below the data by a factor of 
30 as mentioned above. 

It is not surprizing that a good fit of the data can be obtained with two con- 
tributions to the fit that scale as 1/pf and 1/pf at large p t . It is not at all clear 
how accurate the fitted matrix elements are, however, as any effect that modifies 
the slope of the p t -spectrum would affect the relative weight assigned to the two 
matrix elements (Of ( 3 Si)) and Mf^S^^P^). Some uncertainties that enter 
the theoretical prediction can be enumerated: 

1. One may vary the parameters of the calculation such as a s , the renor- 
malization and/or factorization scale, as well as the parton distribution set. A 
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Figure 5: Fit of color octet contributions to direct J/ip production data from 
CDF (-y/i = 1.8 TeV, pseudo-rapidity cut \rj\ < 0.6). Theory: CTEQ4L parton 
distribution functions, the corresponding A 4 = 235 MeV, factorization scale /i = 
(p 2 + Am 2 ) 1 / 2 and m c = 1.5 GeV. Taken from RefM 
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Figure 6: Same as Figure |5] for prompt ip' production. 
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CTEQ4L GRV (1994) LO MRS(R2) 



m 3 j/ /(^ 8) , 3 Po (8) ; 



(8)n 



1.06±0.14t£ : °jj 
4.38 ± 1.15^ 



0.44 ±0.08: 



-0.43 
-0.24 

1.80 ± 0.56lg;i 



1.12 d= 0.14±g;§§ 



1.40 ± 0.22±J;fi 



3.90±1.14±1$ 10.9 ±2.071?;^ 
0.46 ± 0.08±g;i 0.56 ±0.1lig;g 
1.60 ±0.5llg:S 4.36 ±0.96^ 



Table 4: NRQCD matrix elements in 10~ 2 GeV 3 . First error statistical, second 
error due to variation of scale. Ratio of ip' to J / ip fixed. Parton densities from 
RefE. 

different value in a s modifies the over-all normalization and also the slope of the 
Pi-distribution, because the coupling is evaluated at ft = \jAm 2 c + pf . 

The parametrization of the gluon density affects the prediction in a system- 
atic way. Roughly, we may estimate the typical proton momentum fractions of the 
gluons that participate in the hard collision as x g ~ 2\Jpt + iraf/y^i for a given 
value of pt- A gluon density with steeper small-x behaviour therefore steepens the 
Pi-spectrum for the 3 Si component. As a consequence, the M^ 5 ( 1 Sq 8 \ 3 Pq 8 ' ') com- 
ponent whose magnitude depends crucially on the data being somewhat steeper 
than the 3 S[ 8 ^ component, decreases for a steeper gluon density. 

The combined effect of a lower a s and flatter small- x gluon distribution is 
clearly seen in Tab. [| by comparing the value of M 3 ^ 5 ( 1 ^ 8) , 3 P (8) ) obtained with 
MRS(R2) distributions with the one obtained from CTEQ4L or GRV. In contrast, 
the value of (Og^( 3 Si)) is rather insensitive to the parton distribution set, but 
absorbs most of the uncertainty in overall normalization obtained by varying the 
factorization scale fx within a factor two about \J Ami + Pt ( see the second error 
in Tab. |4j). From these variations alone one can conclude that none of the two 
matrix elements is determined to an accuracy better than a factor of 2 and that 
M3 5 ( 1 Sq 8 \ 3 Pq 8 ^) is particularly sensitive to any effect that influence the slope of 
the ^-distribution. 

2. The variation of renormalization and factorization scale gives some insight 
into the magnitude of higher-order corrections in a s only for those contributions 
that have the same kinematic dependence as the leading-order term. In high- 
energy processes or at not so large transverse momentum this may be rather 



26 



misleading. 

Looking at Fig. |5| or ^ we may wonder why at p t ~ 2m c the contribution from 
color octet 3 Pj and cc pairs exceeds the color singlet contribution by about a 
factor 20, although from the scaling rules the color octet contribution would have 
been expected to be suppressed by v 4 ~ 1/10. Since the ratio of matrix elements 
is indeed Mf^S^ , 3 P^ } )/ (Of ( 3 Si)) ~ 4 • KT 2 , the short-distance coefficient of 
the octet term must be three orders of magnitudes larger than the coefficient of 
the singlet term. Examining the short- distance coefficients at p t = 2m c , one finds 
approximately 

da/drips?] 81 (2^r' (35) 

where s is the partonic cms energy. The constant 81 is a large numerical factor, 
partly accidental, partly related to color factors and the larger number of P-wave 
intermediate states. The parametric ratio s 2 /(2m c ) 4 follows from the fact that 
diagrams with t-channel gluon exchange do not contribute to the leading order 
color singlet amplitude. Thus all propagators in the color singlet amplitude are 
off-shell by s, while the octet amplitude is dominated by t-channel exchange that 
allows propagators to be off-shell by only 4m 2 . Because the gluon density favors 
small- x gluons, s is not tremendously large at the Tevatron. Nevertheless, at 
p t = 2m c , the lower kinematic limit specifies s > (l + v^2) 2 (2m c ) 2 ~ 5.8 (2m c ) 2 and 
the ratio s 2 /(2m c ) 4 provides a further significant enhancement of the amplitude. 
This rough estimate approximately coincides with the factor 10 3 estimated from 
Fig.0 

It follows from this discussion that the next-to-leading a 4 color singlet con- 
tribution to which t-channel exchange diagrams contribute should be strongly 
enhanced as well and lead to large i^-factors that increase with increasing trans- 
verse momentum, rather similar as in photoproduction.il At p t ^> 2m c , this 
contribution falls as l/p^, slower as the lowest order contribution (but still faster 
than the fragmentation component (|30D ) and rather similar in shape to the 3 Pj 
and octet components. The value of Mf^S^^PQ 8 ^) will therefore decrease 
to compensate for the additional singlet contribution. 

Another source of potentially large higher-order corrections is related to initial 
state radiation. These, as does intrinsic transverse momentum, would primarily 
modify the ^-distribution in the small-p t region, up to at least several GeV. Ini- 
tial state radiation could be treated either by analytic resummation or estimated 
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by Monte Carlo event generators. A first step in this direction has been taken.0 
The bulk effect seems to be an enhanced short-distance coefficient that leads to a 
decrease of both fitted color octet matrix elements. For a comparison of (0$ ( 3 Si)) 
in Ref.0 and those quoted here and in RefS one has to take into account that the 
evolution of the fragmentation function in the 3 S[ S ^ channel has not been incor- 
porated in Ref.S. Evolution depletes the fragmentation function at large z and 
consequently increases the fitted value of (Of ( 3 Si)) by about a factor of 2. 

3. Apart from higher-order radiative corrections, higher-order corrections in 
v 2 , the parameter of the non-relativistic expansion, can be important close to 
boundaries of partonic thresholds, as discussed in Sect. [27|. Such a situation 
occurs for the color octet fragmentation function that enters the most important 
fragmentation process ([H]) at large p t . At the scale 2m c it is given byi 

D g ^(z, 2m c ) = 6(1 -z) (OSCSi)). (36) 

The delta-function clearly neglects that a fraction Sz ~ v 2 of the gluon momentum 
is transferred to light hadrons in the 'hadronization' of the color octet intermediate 
state. Because of the steep pt-distribution, negligence of this non-perturbative 
softening of the fragmentation function introduces a systematic biasiBO towards 
a too low fitted value of (Cg ( 3 <Si)). 

The momentum taken by the light degrees of freedom (light hadrons) is in- 
corporated into the velocity expansion of NRQCD through operators with cms 
derivatives. Resumming the leading contribution results in0 

D g ^(z,2m c ) = na ^ c) I dy + S(l - z - y + ) /^[ 3 sf ] (y + ), (37) 



c 



where 



£<°lxVr°#(A) + X) 5(y + - [A ■ (iD)] + /k + ) (0(A) + X|^VT o x|0). 



D denotes a (cms) derivative in the quarkonium rest frame, A the Lorentz boost 
that transforms it into a moving frame, '+' the light-cone 'plus'-component of a 
four-vector and k is the four-momentum of the 'hadronizing' color octet cc pair. 
The distribution U[ 3 S[ 8) ](y + ) can be interpreted distribution of light-cone 
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momentum fraction y + taken by light hadrons in the pro cess :fffj 

cc[ 3 S[ s) ] -f tp + light hadrons. (39) 

After resummation the octet fragmentation function is expressed in terms of a 
new unknown function and all predictivity seems to be lost. However, we know 
that /</,[ 3 5i ](y+) should have support mainly in a region y + ~ v 2 . Moreover, it is 
a universal function that could in principle be extracted from other processes not 
necessarily related to fragmentation. Once such a universal function is isolated, 
one can also try to model it and check the consistency of the ansatz in different 
production processes. 

Summarizing this discussion, I think that regarding the numerical values of the 
matrix elements in ([34]) as accurate within a factor of 2 or less for Mf^S^^P^) 
would not be overly conservative. It is worth noting that the color evaporation 
model mentioned earlier can reproduce the Tevatron data rather reasonably M 
This is no surprize for the shape of the ^-distribution, as the CEM includes a 
fragmentation component. It is less obvious that this works (again within fac- 
tors of 2) with the same values of normalization factors required to describe 
fixed-target data. From the point of view of NRQCD such agreement would be 
considered coincidental, perhaps related to the fact that the octet matrix elements 
relevant to both processes all scale equally as v A . 

Prompt S'-wave bottomonium production has also been measured by CDF, 
although a Xft-component could not yet be separated. A comparison with pre- 
dictions can be found in Ref.il. As v A ~ 1/100 for bottomonium color octet 
contributions are less relevant than for charmonium in the region p t < 15 GeV 
of the Tevatron data. Given the uncertainties in the color singlet contribution 
discussed above and the octet matrix elements for bottomonium, the necessity 
of color octet contributions is not yet as firmly established as for charmonium 
production. A realistic description of the p t -spectrum requires the resummation 
of soft gluon radiation, which has not yet been undertaken. 



R While writing this review, I learnt of related ideas by M. Mangano, who introduces a sim- 
ilar distribution function on phenomenological grounds outside the context of the NRQCD 
expansion. 
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3.2 Polarization 



Perhaps the most decisive test of the NRQCD picture of quarkonium production 
will come from a polarization measurement of direct J/ip and ip' at the Tevatron 
at large quarkonium transverse momentum. Recall that at large p t , the produc- 
tion cross section is dominated by gluon fragmentation into ccpS^ ]. It was first 
noted by Cho and WiseS that the cc pair is transversely polarized, because the 
coupling of longitudinal gluons is suppressed by 4m 2 /pf. Furthermore, to leading 
order in the velocity expansion the subsequent transition (p9|) takes place via a 
double electric dipole transition, which does not flip the heavy quark spin. Conse- 
quently, at large transverse momentum, one should observe transversely polarized 
quarkonia. The polarization can be observed in the angular distribution of decay 
leptons from ip — > which can be written as 

ocl + acos 2 6», (40) 



dcosO 



with 9 the angle between the lepton three-momentum in the ip rest frame and the 
polarization axis, the ip direction in the hadronic cms frame (lab frame). Pure 
transverse polarization would imply a — 1. 

Corrections to the asymptotic limit come from three sources, corresponding 
to the three small parameters v 2 , a s /iT and 4m 2 /p 2 that characterize quarkonium 
production at the Tevatron.Q Since spin symmetry is violated by higher order 
interaction terms in the NRQCD Lagrangian, longitudinal polarization can arise if 
the transition (|39]) proceeds via two magnetic dipole transitions. These corrections 
scale as v 4 and do not vanish in the limit of large transverse momentum. According 
to the power counting estimate^ they could reduce a to 0.85. 

The second source of corrections comes from a s corrections to the fragmenta- 
tion function, since radiation of a hard gluon can change the cc pair's polariza- 
tion. This correction decreases with p t as l/(lnp t ). All relevant fragmentation 
functions into longitudinally polarized quarkonia have been computed^ and were 
found to yield about 5% longitudinal polarization at p t = 20 GeV. The correc- 
tions turned out to be rather small, mainly because the fragmentation functions 
can only be softer than the leading term (|36| ) that produces transverse charmo- 
nia. After convolution with the short- distance cross section this leads to sup- 
pression of the higher-order contribution. If added to the first correction, the 

*I am ignoring higher- twist corrections which are suppressed in the small ratio A/pj. 
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Figure 7: a as a function of pt- Taken from Ref. ED 

second still implies that a > 0.65 at p t = 20 GeV. This estimate is based on 
(Ot{ 3 P ))/m 2 c « (Oti^)), consistent with Tab. | (CTEQ4L). 

Below p t ~ 15 GeV the dominant source of depolarization is due to non- 
fragmentation contributions,^ because they have large short-distance coefficients 
as seen above. The size of these depends crucially on the magnitude of the pa- 
rameter M^' 5 ( 1 5'o 8 ' ) , 3 Po 8 ' ) )- With the value taken from Tab. |], the prediction for 
the parameter a in the polar angle distribution is shown in Fig. |7]. The band is 
based on the non-fragmentation corrections and a s corrections to the fragmenta- 
tion functions, but neglects the depolarization due to spin symmetry breaking of 
order 8a ~ 0.1. It is important that the band reflects only the uncertainties in 
the extraction of NRQCD matrix elements due to statistical errors in the data of 
the unpolarized cross section. It does not include any of the theoretical uncertain- 
ties that could systematically affect the extraction of these matrix elements. In 
particular, if Mf 5 ( l S^\ 3 PQ 8 ^) is smaller than assumed here, which is not unlikely, 
the transverse polarization fraction would increase. Nevertheless, at pt ~ 5 GeV, 
Fig. |7| shows that no remnant of transverse polarization may be expected. As 
p t increases, the angular distribution becomes rapidly more anisotropic. If above 
Pt ~ 20 GeV a substantial fraction of transverse polarization is not observed in 
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future (run II) Tevatron data, we will definitely know that NRQCD, or at least 
spin symmetry, is not applicable at the charm mass scale. If it is observed, the 
color evaporation model will have lost most of its remaining appeal. 

4 Quarkonium production at fixed target 

It has long been known, and highlighted most recently in Refs. , that the color 
singlet model does not give a satisfactory description of all quarkonium production 
cross sections at fixed target energies. The discrepancies arise for direct J/ip and 
Xi, both of which can not be produced by gluon-gluon fusion at leading order in 
the color singlet channel. When the over-all normalization of the color singlet cross 
section is adjusted to the data, one also predicts a significant fraction of transverse 
polarization, a ~ 0.25 in the notation of fllPp and in the Gottfried- Jackson frame. 
No polarization is observed. As compared to the Tevatron '^'-anomaly' these 
discrepancies are less dramatic and the theoretical uncertainties could always be 
stretched, and higher-twist effects invoked, to make them appear even less incisive. 
Nevertheless, these discrepancies must be reassessed in the light of the NRQCD 
factorization approach. 

4.1 Cross sections 

Fixed target quarkonium cross sections have been calculated by several collabo- 
rations.0Hi§ The spin and color states that can be produced in parton-parton 
collisions at order ot 2 s are as follows: 

gg - cc[n): n = l S^\ *P§f , l P{ 8 \ D-waves, . . . 

qq — > cc[n] : n = 3 S[ 8 \ D-waves, .... (41) 

At order a 3 , after hard gluon radiation, no stringent helicity and color constraints 
remain. Adopting a value for (Of ( 3 Si)) on the order of those shown in Tab. |], the 
quark annihilation channel turns out to be insignificant at energies far enough from 
the cc threshold such that X1X2 ~ 4m^/s <C 1 for the product of parton momentum 
fractions. The cross sections are then dominated by gluon-gluon fusion for both 
pion and proton beams. Keeping the intermediate cc states that are leading 
according to the scaling rules of Tab. ||, the direct J /if) and ip' production cross 
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Figure 8: Total (solid) and singlet only (dotted) ip' production cross section (xp > 
0) in proton-nucleon collisions. The solid line is obtained with Mf ^S^^Pq 8 ^) = 
5.2 • l(r 3 GeV 3 . 

section is given by 



&{S9 -> i>') 

5ii 3 a 2 
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where z = (2m c ) 2 / {SX1X2), is the center-of-mass energy and a s is normalized 
at the scale 2m c . The third line gives the color singlet cross section that enters 
at order a 3 . Using spin symmetry to relate (Ot( 3 Pj)} = (2 J + 1) (C^( 3 P )), the 
leading color octet contributions are proportional to Mf^sjpfPjP) defined as m 



A comparison of the predicted energy dependence and normalization of the 
cross section (including the qq annihilation channel) with data is shown in Fig. |8| 
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Figure 9: J/ip production cross sections (xp > 0) in proton-nucleon collisions. 
Solid line: Fit to the total cross section including radiative feed-down from the 
Xcj and if)' with '{^ ?Pf>) = 3.0 • l(r 2 GeV 3 . Dashed line: Direct J/tf; 
cross section. Dotted line: Direct J/if> production, color singlet only. 

and P for if)' and J / ip production, respectively. The solid curves are obtained with 

MfCsjPfP®) = (Of'CSo)) + — 2 (Ot'( 3 P )) = 0.56 ■ 10- 2 GeV 3 

M 7 J/ ^( 1 4 8) , 3 P (8) ) ee (Of'CSo)) + (O J 8 ^( 3 P )) = 3.0 ■ 10- 2 GeV 3 , (43) 

with m c = 1.5 GeV, which is assumed in this analysis. The color singlet contribu- 
tions to the direct production cross section is shown separately. Note that since 
the qq annihilation channel is insignificant, the P-wave feed-down incorporated 
for J /if> is entirely color singlet, because relativistic corrections to xi production 
have been neglected. I will discuss % production in more detail in Sect. |4.2| . 

Turning first to if)' production, one notes that given the uncertainties inherent 
to leading order QCD predictions (no i^-factors have been introduced), their 



34 



renormalization scale-dependence, as well as the dependence on m c , the relevance 
of color octet contributions is far from obvious. This is mainly due to the fact 
that the color singlet cross section is dramatically enhanced compared to Ref.,@ 
because it is expressed in terms of 2m c = 3 GeV rather than the significantly larger 
ip' mass. Expressing the cross section in terms of quark masses implies a product 
X1X2 smaller by 50% than its physical value and enhances the cross section as a 
consequence of the rising gluon density at small x. As noted earlier, this intrinsic 
ambiguity of a leading order calculation in v 2 would be alleviated by resumming 
higher-dimension operators. In the present case it is likely that the effect would 
be a reduction of the color singlet cross section, indicating more stringently a 
missing additional production mechanism. The short-distance coefficients of color 
octet production would be likewise affected by implementing the physical phase 
space constraints. In fact, as the invariant mass of the final state containing 
and light hadrons is even larger than m^r, an even more significant reduction 
would be expected. A fit of the normalization to data would then require a larger 
Mf(iS®*pW) than the one quoted in ( [13]) , which would be desirable in view of 
the apparent discrepancy between ( pE3] ) and (|3H|). Given these uncertainties, we 
can conclude that the sizes of color octet contributions as deduced from Tevatron 
and fixed target data are consistent with each other as far as order of magnitude 
is concerned and that one could not have hoped for more than that. 

For J/ip, Fig- we have a clear deficit in the color singlet direct production 
cross section. The difference to if)' comes from the fact that mj/^ is not much 
larger than 2m c . While this does not remove the ambiguity with the choice of 
m c , one would at least seem to be already closer to the physical phase space. 
If we assume that the P-wave feed-down is determined by leading order color 
singlet cross sections, an absolute prediction is obtained for the direct J/ip cross 
section, once the unknown parameter M./^ ^S^^Pq 8 ^) is fitted to the total J/ip 
cross section. One finds that a{J /VOdW 'cr(J Jijj) ~ 0.6 in good agreement with 
experiment. On the other hand this ratio comes out a factor 3 too small with color 
singlet contributions only. The result <r{J/ip)&iJ 'cr(J '/ip) ~ 0.6 is still somewhat 
problematic, because it relies on the assumption that the x feed-down is accurately 
reproduced and understood. But, as discussed in the following subsection, the 
Xi/Xi composition of the x feed-down assumed here conflicts with experiment 
in that there are too few xi predicted. Thus there should be an additional, but 
neglected x production mechanism that would then upset the nice agreement 



35 



found for a(J /ip)dh/ 'cr(J /ip) (and make the color singlet model look even worse). 
This, however, can be compensated by taking into account the larger mass of \ 
states in the same way as for if)'. Thus, it seems that color octet contributions 
in fixed target production are natural candidates to remove large discrepancies in 
over-all normalizations, but the remaining uncertainties preclude a straightforward 
test of the universality of color octet matrix elements by comparison with Tevatron 
data. In particular, one can not exclude the possibility of significant 'higher-twist' 
contributions that would have to be added on top of color octet contributions. The 
numerical situation will be considerably improved by the soon available next-to- 
leading QCD calculation of all relevant octet and singlet production channels. For 
P-waves this analysis has been reported in Ref.0 

A comparison of cross sections in experiments with pion beams, using ( pE3"D 
extracted from proton beam experiments shows that the former are then system- 
atically under-predicted by a factor of about 2. The resolution of this discrepancy, 
already known in the context of color singlet model predictions,^ is clearly outside 
the scope of NRQCD, which assumes factorization of the hadronic initial from the 
hadronic final state. 

The analysis of bottomonium S-wave cross sections suffers from the lack of 
sufficient experimental data, and in particular the absence of any experimental 
information on Xb cross sections. Using color singlet wave functions and an esti- 
mate of color octet matrix elements, one finds0 a larger indirect contribution to 
S- wave cross sections than for charmonium. The observation of T(3S) comparable 
to T(ljS') in the E772 experiment^ then indicates the existence of yet undetected 
Xb(3P) states below the open bottom threshold. 

In addition to integrated cross sections, rcpSand p$ distributions may be 
analysed. Since the xf distribution is determined by the gluon flux, the distri- 
butions are identical for color singlet and octet contributions and do not lead to 
further insight. Note that at large xp the NRQCD expansion breaks down and 
higher-twist corrections become large.0 A realistic comparison of p t distributions 
in the range of available fixed target data requires yet accounting of soft gluon 
effects. 
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Experiment 


beam/ target 


y/s/GeV 


R 


E673 


pBe 


19.4/21.7 


0.24 ± 0.28 


E705 


pLi 


23.7 


0.09 


E771 


pSi 


38.8 


0.34 ± 0.16 


WA11 


7rBe 


18.63 


0.7 ± 0.2 


HvO ( 6 


i ^ 

TYDC 


lo.oo 


U.yo ± U.o4 


E673 


7rBe 


20.55 


0.9 ±0.4 


E705 


7rLi 


23.72 


0.53 


E672/706 


7rBe 


31.08 


0.57 ±0.19 



Table 5: Experimental data on the X1/X2 ratio R. Data compiled from Refs.EJ. 
In view of the experimental errors no attempt has been made to rescale older 
measurements to account for the latest \ branching fractions. 

4.2 The Xi/%2 ratio 

As noted in (|4~1~D, the leading-order gluon-gluon fusion process can not produce a 
spin-1 P-wave state. Consequently, at energies where qq annihilation is irrelevant, 

R = ^4 = o( — ) ~ 0.05 (color singlet), (44) 

with the numerical value from Ref. 0. The datai compiled in Tab. || clearly show 
more copious Xi production. 

In the context of the color singlet model, the Xi production cross section has 
sometimes been boosted by taking into account an exclusive quark-antiquark fu- 
sion mechanismil qq — > Xcj at order af, see e.g. the cross sections displayed in 
Fig. 16 of RefS. As this mechanism produces Xcj states in the ratio 0:4:1, it 
raises R at low energies, especially in tcN and pN collisions. The justification for 
keeping this higher-order contribution derives from the presence of two infrared 
logarithms, so that this contribution is of order c^ln 2 (m c /A), where ln(m c /A) is 
interpreted as a non-perturbative parameter. This bears some similarity with the 
infrared logarithm in qq Xcj + 9 which, as discussed, could be interpreted as 
qq — > cc[ 3 S'j 8 ' ) ] — > XcJ + A in the framework of NRQCD factorization. This similar- 
ity is superficial, though, because the double logarithm above is not regularized by 
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the quarkonium binding energy and can not be absorbed into a NRQCD matrix 
element. Thus, if it existed, it would contradict NRQCD factorization. However, 
it seems clear that the origin of the double logarithm is the exclusive final state 
considered. The relevant cut diagram at order a 4 contains other cuts related to 
interfering qqg final states, which would cancel the infrared logarithms in agree- 
ment with NRQCD factorization. Thus, there is no infrared enhancement in the 
quark-antiquark fusion mechanism and it can not help to raise the X1/X2 ratio. 

The next candidates are relativistic corrections in the velocity expansion.BH 
Indeed, since they dominate direct J/ip production, that suffers from the same 
suppression as xi i n t ne color singlet gluon-gluon fusion channel, they should 
naturally also be large for xi production. The contributions that scale as v A 
relative to the leading order contribution for P-wave production are listed in 
Tab. |3|. (In each channel, there exist v 2 corrections due to operators with spatial 
derivatives squared.) To get an estimate with no pretence of accuracy, let me 
consider only the octet 3 Pj> intermediate states. They yield 

q( Xl )from 3 Pf 15 (Of( 3 P 2 )) + 15/4(Of ( 3 P )) 1 
( T(x2)from3p« 8 <0f( 3 P 2 )> ~ 3' 1 J 

where I assumed that (O^ 1 ( 3 P 2 )) / (O^ 2 ( 3 P 2 )) ~ 1/10 according to the scaling v 4 ~ 
1/10. If we multiply this ratio by 2 to account for the intermediate S- and P-wave 
states, we obtain a xi production cross section an order of magnitude larger than 
([44]). Since the new production channels contribute also to X2 production, let me 
assume that they contribute to Xi an d Xi i n the ratio 3:5. Then 

R « 0.3. (46) 

I emphasize again that this is a very crude estimate.[j] The estimate is in agreement 
with the X1/X2 ratio measured in pN collisions, but is lower than the measurement 
in ttN collisions. If the trend of larger R for pion collisions exhibited by the 
measurements shown in Tab. |^ is confirmed, the real problem would no longer 
be the X1/X2 ratio per se, but the difference between proton and pion beam 

tThe short-distance coefficients of all v 4 suppressed intermediate states have recently been cal- 
culatedS Nevertheless, the significant number of new matrix element makes it difficult to 
improve on the crude estimate given here. The estimate of R given in Ref.0 is larger than the 
one quoted here, because for unknown reasons the authors choose to omit the dominant color 
singlet channel in \2 production. Accounting for it makes their estimate consistent with (|4q). 
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experiments. Unless there is an unknown enhancement of the qq annihilation 
channel, such a difference would violate the initial/final state factorization of 
NRQCD and would have to be attributed to a 'higher-twist' effect. Note that as 
long as the qq annihilation channel can be neglected, the X1/X2 ratio should be 
approximately energy- independent. 



4.3 Polarization 

Polarization measurements have been performed for both if) and if)' production 
in pion scattering fixed target experiments.ilS Both experiments observe an 
essentially flat angular distribution in the decay if) — > {if) = J/if),if)'), 

oc 1 + a cos 2 6 , (47) 



dcos9 



where the angle 9 is defined as the angle between the three-momentum vector of 
the positively charged muon and the beam axis in the rest frame of the quarko- 
nium. The observed values for a are 0.02±0.14 for if)', measured at y/s = 21.8 GeV 
in the region xp > 0.25 and 0.028 ± 0.004 for J /if) measured at yfs = 15.3 GeV in 
the region xp > 0. In both cases the errors are statistical only.fj] 

Compared to these measurements, the color singlet contributions alone yield 
a m 0.25 for the direct S- wave production cross section.@ The polarization yield of 
octet contributions has been considered in RefR|§] The ccpS^] intermediate state 
yields pure transverse polarization, but is insignificant for the total cross section, 
because it occurs only in qq annihilation. The cc[ 1 S u 8 ' ) ] intermediate state results 
in unpolarized quarkonia while the cc[ 3 Pj 8 ' ) ] intermediate states prefer transverse 
to longitudinal polarization by a ratio 6:1. The longitudinal polarization fraction 
is therefore proportional to Mf ( 1 Sq 8 ' ) , 3 Pq 8 ' ) ), different from the combination that 
enters the cross section summed over all polarizations. Constraining (Og( l So)) 
and (Of ( 3 -Po)) to be positive, the range 

0.15 < a < 0.44 (48) 

Contemplating Fig. 11 of Ref.!^ I am strongly tempted to assume that the error is misprinted 
and that the measurement should read 0.028 ± 0.04. 

S Polarization has also been treated in RefR However, in the first reference of Refill and the 
preprint version of the second it was assumed that the ^g 8 ' production channel yields pure 
transverse polarization rather than no polarization, which leads to results significantly different 
from RefR This has been corrected in the published version of the second reference. 
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is obtained for if}' production at ^/s = 21.8 GeV. (The energy dependence is rather 
mild.) The lower bound is attained if (Of ( 3 Pq)) = and could be as low as 



0.08, since, as discussed in Sect. |4.1| , the color singlet fraction in ip' production 
is likely to be lower than that shown in Fig. [8|. Thus, the prediction can be con- 
sistent with data within experimental errors, if (Of ( 3 Po)} is small compared to 
(Of (^o)). Such a scenario would be realized in case that the multipole suppres- 
sion of (O'f (^o)) is not effective (A ~ 1 in Tab. 0). 

The discussion for direct J /if) production follows the ip' case. To compare 
with the measurement, the polarization yield from x an d ip' feed-down has to be 
computed. One then obtains 

0.31 <a< 0.62, (49) 

the larger value of a being related to the fact that only color singlet mechanisms 
to x production were considered (so that the feed-down is almost purely X2, m 
conflict with experiment) and that \2 is produced only in a helicity ±2 state in 
gluon-gluon fusion at leading order. The feed-down therefore yields an almost 
purely transversely polarized component to the J /if) cross section. Eq. (^) is 
incompatible with the observed flat angular distribution. As a prediction it should 
however be taken with due caution. Taking the estimate of Sect. [4.2| literally, only 
one third of the x feed-down is due to X2 produced in a color singlet state, while the 
polarization yield of the remainder is unknown. Thus, it does not seem excluded 
that the total x feed-down is largely unpolarized, in which case a can be as small 
as 0.1. To resolve this issue, a measurement of x polarization would be desirable. 
In particular, one would like to know to what extent Xc2 is produced with helicity 
±2. Meanwhile, the situation remains unclear and one may speculate on the role 
of higher-twist final state interactions without impunity. 

The comparison with polarization data, as well as the compatibility of matrix 
elements determined from large-p t and fixed target production, Eqs. ([34]) and ([43]) 
could be ameliorated, if, as has been suggested,^ (Of ( 3 Po)) is actually negative. 
This is indeed possible in principle for the renormalized matrix element, because 
of operator mixing, and in particular a quadratic divergence that would appear 
in schemes with an explicit factorization scale. However, phenomenologically, a 
negative value is unacceptable. The quadratic divergence implies that in any 
process (Of ( 3 Po)) appears as a combination 

(OfCPo^ + kasfi^Of^S,)), (50) 
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where \i denotes the factorization scale and the constant k follows from separating 
and including only the soft gluon contribution to the color singlet cross section. 
As this combination enters a physical cross section, it must be positive in any 
factorization scheme. On the other hand, in all phenomenological applications, 
the higher order color singlet piece is not included into the theoretical predic- 
tion or is argued to be smaller than the octet contribution. Consequently, the 
phenomenologically extracted value of (Of ( 3 Po)) alone must satisfy the positivity 
constraint. If a fit prefers this matrix element to be negative, it must be attributed 
to uncertainties in either theory or data. 

4.4 Beyond NRQCD 

As emphasized in Sect. |2.2| , the NRQCD factorization approach is leading twist 
and neglects potential interactions of the intermediate heavy quark pair with the 
beam or target remnants, as well as any production mechanism that involves multi- 
parton interactions in the initial state. Quarkonium production cross sections in 
fixed target experiments, where most quarkonia escape at small angle with the 
beam axis, would be especially vulnerable to such corrections. In this subsection, 
I list some of such phenomena, that can not be described by NRQCD. For more 
comprehensive discussions and references, see the overviews given by BrodskyS 
and by Hoyer.il 

Nuclear dependence. Charmonium production on nuclear targets shows sup- 
pression with increasing nucleon number, which is parametrized by o = <JoA a , 
where a ~ 0.9 — 0.92 for cross sections integrated over positive xp. Since a ~ 1 
for open charm production, this nuclear dependence can not be attributed to 
shadowing in parton distributions. NRQCD can not account for such a nuclear 
dependence due to rescattering of a cc pair with the nucleus. Since for bottomonia 
a ~ 0.97, nuclear suppression is consistent with being higher twist. Note that in 
the nucleus rest frame, the charmonium formation time is long in a high energy 
collision and charmonium forms only after the cc pair has left the nucleus. Since 
most J/ip are now believed to originate from a cc pair originally produced in a 
color octet state, the rescattering cross section is that of an octet rather a singlet 
cc pair. A further consequence is that despite their different sizes, the nuclear 
dependence should be the same for J/ip and if)' as is indeed observed. 
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Comover interactions. A cc pair that moves with nearly equal velocity as light 
quarks from the hadron remnants may find it preferable to produce open charm 
rather than charmonium. Comover interactions responsible for this effect would 
be expected to be largest in the nuclear fragmentation region (xp < 0), because 
of the larger number of potential comovers. Consequently, a should decrease for 
xf < 0, an effect that is observed for both charmonium and bottomonium. 

Large-XF phenomena. The region xp — ► 1 is particularly interesting, because 
the usual hard scattering picture breaks down in this region.0 One way to vi- 
sualize this is to note that all partons of the projectile have to be correlated in 
order that all momentum can be transferred to a single parton, that subsequently 
transfers all momentum of the projectile to the cc pair. Such a correlation is 
rather short-lived and at large xp its lifetime can approach the hard interaction 
time scale set by l/m c . A hard collision at large xp then takes a snapshot of a 
compact Fock sate component of the projectile, for which multi-parton interac- 
tions are not suppressed. Such effects may explain the turn-over to longitudinal 
polarization at large xp seen by the Chicago- Iowa-Princeton collaboration, and 
also a decrease of a towards 0.7 at large xp, since the soft parton in the projec- 
tile stopped after transferring its momentum to the cc might then scatter only 
from the surface of the target nucleus. If the hadron wave function contains an 
intrinsic charm component (which it certainly does at some level), it would also 
naturally show itself at large the intrinsic cc pair preferentially carries most 

of the projectile momentum. Charmonium production due to intrinsic charm is 
suppressed as A 2 /m 2 (but can be sizeable in narrow kinematic regions) and is not 
included in NRQCD. 

5 Conclusion 

In this lecture I have tried to convey a topical overview on non-relativistic QCD as 
an effective theory and some of its applications. As any effective theory, NRQCD 
relies on the presence of different mass scales, here the quark mass tuq and the 
quark's typical velocity rngv and energy niQV 2 in a non-relativistic bound state. 
Assuming perturbation theory at the scale tuq, NRQCD allows us to organize 
a calculation as an expansion in a s (mq) and v 2 . NRQCD shares many features 
with heavy quark effective theory and yet has a much richer structure owing to 
the existence of several low energy scales tuqv, mQV 2 and A. As a consequence the 
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power counting in the effective theory is also complicated and one cannot rely on 
the dimension of operators alone. The velocity scaling rules of NRQCD are simple 
and unique in the Coulombic limit mqv 2 ^> A. If this limit were realized, our use 
of NRQCD would mimic non-relativistic QED and the bound state properties of 
quarkonia would be amenable to a perturbative treatment. As A is close to any 
of the two other bound state scales for charmonium, the velocity scaling is more 
delicate (and perhaps the most debated issue in NRQCD) and must in principle 
be assessed for every quarkonium state anew. 

To describe quarkonium production, NRQCD has to be amalgamated with the 
concept of factorization in perturbative QCD. This opens tremendous perspectives 
on phenomenology, of which only a fraction could be reported here. New insights 
into quarkonium production followed from recognizing the importance of color 
octet and fragmentation production mechanisms, both of which are crucial in 
accounting for the charmonium cross sections observed in the Fermilab Tevatron 
experiments. 

Nevertheless, NRQCD is not a 'Theory of Everything' in quarkonium produc- 
tion. The factorization formula for quarkonium production can be compared to 
the factorization theorem for Drell-Yan processes, with higher twist corrections 
due to, for example, initial and final state multi-parton interactions neglected 
in both cases. Such corrections may be more or less severe dependent on the 
production process, but would affect fixed target and forward photoproduction 
in particular. The comparison of leading twist predictions with fixed target data 
showed significant improvement in all aspects compared to the color singlet model, 
although some features such as x production and ip polarization could be made 
just consistent with data only upon straining the theory to its limits, which does 
not appear entirely satisfactory. In any case, even without higher twist correc- 
tions, theoretical uncertainties remain appreciable and complicate a verification of 
the universality of long-distance matrix elements in NRQCD. It seems, at least in 
my opinion, that at this moment there is no discrepancy in any process, including 
those not discussed here, that could not be attributed to one or another difficulty 
in the theoretical prediction, and which thus is understood, even if it is not easily 
remedied. As concerns fixed target experiments, much clarity could be obtained 
from measurements of x polarization, polarization of T(nS) and separate Xb cross 
sections. 

For any theory, the question 'Can it be wrong?' is at least as important as 
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'Is it correct?'. Since, in the asymptotic limit itlq — > oo, NRQCD is evidently as 
correct as QCD itself, one should ask whether it could generally be inapplicable 
at the charmonium scale. Here, at last, the situation is clear: ip' or direct J/ip 
polarization at the Tevatron at large transverse momentum will evidence dramatic 
success or failure of NRQCD for charmonium production. 
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